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Abstract 

We describe an experimentally straightforward method for preparing an entangled W state of 
at least 100 qubits. Our repeat-until-success protocol relies on detection of single photons from 
collective spontaneous emission in free space. Our method allows entanglement preparation in 
a wide range of qubit implementations that lack entangling qubit-qubit interactions. We give 
detailed numerical examples for entanglement of neutral atoms in optical lattices and of nitrogen- 
vacancy centres in diamond. The simplicity of our method should enable preparation of mesoscopic 
entangled states in a number of physical systems in the near future. 
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The preparation of entangled states is a fundamental task in quantum information 
processing (QIP). A wide variety of physical systems exhibit long-lived quantum coherence 
and are potentially useful for QIP, but lackj^ 
that permit entangled state generation 



he strong, controllable qubit-qubit interactions 
0,0. Various measurement-based quantum 



computing schemes have demonstrated that the repeated preparation of specific entangled 
states is an adequate substitute for unitary entangling operations 0, |j, 5|. Here we 
propose a method for preparing the entangled W state of at least 100 qubits without 
any direct or indirect qubit-qubit interactions. The W state of N qubits can be written 
\W N ) = 1 100. . .0) + |010. . .0) + |001. . .0) + . . . + |000. . . 1). This class of states exhibits 
unusually robust nonclassical properties p, Q] and can be used for quantum teleportation 
and secure communication j^, 9]. Our method is applicable to a wide variety of qubit 
implementations in which the qubit states exhibit spontaneous emission. Detection of a 
single spontaneous emission photon from a coherently excited assembly of qubits indicates 
successful state preparation. The highly directional character of the collective emission 
in free space allows for a robust and efficient detection setup. The simple experimental 
apparatus required for our method should enable preparation of mesoscopic entangled 
states in a number of physical systems in the near future. 

Our method allows straightforward entanglement generation for many attractive QIP 
candidates that are currently experimentally challenging. Entanglement generation usually 
requires engineering indirect qubit-qubit interactions through an auxiliary quantum degree 
of freedom common to all the qubits, a so-called "quantum bus". The striking success of 
ion-trap QIP relies on the ion-ion Coulomb repulsion as a naturally occurring quantum bus 
[lol |. In most cases, though, the quantum bus itself requires considerable effort to construct, 



and the usual technique is to invoke the highly challenging technology of strong- 



coupling 



cavity QED, as in proposals for QIP with neutral atoms [ll], with quantum dots 12], and 



with NV-centers 13| . In contrast, our reliance on detection of collective emission ensures 



the coherence of the entangled state, even for emission in free space. 



D 



Collective emission processes have recently found use in generation of single photons 

ement between collective excitations of two 



l4 . 1 1 51 ] and in the preparation of entang 



separate ensembles of atomic qubits [la, Il7j . The repeat-until-success protocols used in 
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those works are the inspiration for our method, which generates multipartite entanglement 
of the mesoscopic number of individual qubits making up the ensemble, rather than using 
a collective excitation as a single qubit. The multipartite entanglement available from our 
method is more powerful than entanglement of two collective-excitation qubits, and simpler 
to implement than entanglement of large numbers of collective excitations. Our method 
also avoids the requirement on collective-excitation protocols to use optically dense qubit 
ensembles, which can restrict applications to inconveniently high numbers or densities of 
qubits. 



Collective spontaneous emission 

Consider an assembly of N two-level atoms with transition wavelength A and decay time 
t. The atoms are contained in a region with dimensions much smaller than cr, but much 
larger than A. The atoms are initially prepared in the ground state and interact with a series 
of resonant plane- wave laser pulses of wavevector k^. The atomic dynamics is conveniently 



described using the collective operators 
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parametrized by wavevectors k, where the are the Pauli operators of the jth atom. 
In this formalism, a laser pulse with area 7r/2 and phase 0^ acts as the operator 
U n /2{4>L) = exp[z(7r/4)e~^ L J*± +h.c.]. The radiative dynamics is generally fast compared to 
the motion of the atoms, so we take the atoms to be in eigenstates of the position operator x. 

The radiation produced by atomic decay is described by the atom-field interaction Hamil- 
tonian H mt = H^2? gza?J+(k) + h.c, where is the annihilation operator for a plane-wave 
field along k. The angular dependence of gives rise to the single-atom dipole radiation 
pattern. The intensity distribution of the atomic radiation at short times is found by Fermi's 



Golden Rule to be 



19] 



3 



I(k) = I (k)-(l + C(k)) (2) 
C(k) = -\^M-Kk - k L ) ■ xj)\ 2 (3) 

where Io(k) is the single-atom radiation pattern. The coherent term ((k) is approximately 
equal to N over a narrow range of k ~ fc^, so the radiation pattern contains a narrow peak 
in the forward direction. The power in the forward peak can be a significant fraction of the 
total power. 



Conditional dynamics and entangled-state preparation 

Motivated by this examination, we study the atomic dynamics conditional on emission 
of a single photon along k^. Suppose the atoms are excited with a 7r/2 pulse with phase 
4>l, radiate freely for a time T^et much shorter than the time needed for emission of a single 
photon, and are de-excited with another tc/2 pulse with phase (fti + Tt. If no photon is emitted 
during T det , the state of the atoms is the same before and after the operations. However, if 



a photon is emitted along k L , the pulse sequence implements the operator 
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M det (0 L ) = U W/2 ^ L + tt) J1 l U n/2 (<p L ) (4) 

= l4 L ~ \jl h e~** L + -^J z e- i4>L (5) 
2 2 y2 

where we have assumed that the 7r/2 pulse durations are always much shorter than T& et . 

Successful single-photon detection during the detection period Td e t implies that the atoms 

are in the state M dct (4> L )\gg ...g). 

We show the experimental arrangement for entangled state preparation in Fig. [lj In 
one trial of the repeat-until-success protocol, we detect the atomic radiation along the for- 
ward direction ki during a pair of pulse sequences of the type that implement Md e t- The 
laser phase is shifted by 7r between the two pulse sequences. Before detection, the emitted 
radiation passes through an unbalanced Mach-Zehnder interferometer to induce construc- 
tive interference between the two emission time-bins. Detection of a single photon in the 
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FIG. 1: Experimental arrangement for preparing the W state of N qubits. Two laser pulse se- 
quences generate small amounts of collective spontaneous emission from the assembly of atoms 
(green). An unbalanced Mach-Zehnder interferometer and subsequent temporal gating of the over- 
lapped emission (red) allow selective detection. A click at the detector (red) implements the 
entangling operator M ent . 

interference time-bin implies that the operator 



M ent = ~(M det ((f) L ) + M dct ((f) L + tt)) 



;(4 



Jl L e 



-2icj> L \ 



(6) 
(7) 



has acted on the atomic state. Since annihilates the initial state, we find the final state 
J+\gg . . . g) . This is just the first excited Dicke state along ki for an extended atomic 
ensemble. Regarding each atom as a qubit over the states \g), \e), we see that the final 
state is related by single-qubit phase factors to the W state of N qubits [6], which we write 
|Wjv). As long as the qubits remain motionless on the scale of a wavelength, we can ignore 
the additional phases. 



Sources of error in the entangled state 

As in other repeat-until-success quantum protocols, entanglement purification is built 
into our scheme by the conditional detection [16( . We can analyze the entanglement of the 



final state using the entanglement witness W = (1 — 1/N)t — \Wn)(Wn\ |20|, [21J. For any 
density matrix p, the condition Tr[Wp] = 1 — 1/N — F < implies that p is entangled. 
The infidelity 1 — F is never greater than the sum of the various error probabilities, so the 
final state is certainly entangled as long as the sum of the error probabilities is less than 1/N. 
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For atom number densities less than 1/A 3 , the total emission probability P em i S during 
each trial considerably exceeds the probability of successful detection P^ et . However, 
this inefficiency need not affect the fidelity of state preparation. If a photon is detected, 
indicating a successful trial, there is a probability P C mis that a second photon has also been 
emitted during the detection period. The final state resulting from the double-emission 
event is not the desired entangled state. This error source can be reduced to any desired 
level by shortening the detection time, at the cost of a lower success probability per trial. 

A similar error source arises from the uncertainty in the wavevector of the detected 
photon k det . Mismatch between fc det and k^ induces an additional phase <pj = (k^et — &l) • %j 
on the jth atom, so that we obtain the state \W' N ) = [Hj exp(— i(fij(Jz )]\Wn). Averaging 
the fidelity | (W n\W n ) | 2 = ((k)/N over the cone of detected wavevectors, we find the 
error in the final state to be 1 — (/N. We can reduce this source of error by restricting 
the angular range of detection to keep ((k) ~ N, again at a cost of lower success probability. 

Finally, successful entangled state preparation relies on the initialization of the atomic 
state to \gg . . . g) before each trial. An unsuccessful trial implements the operator 
1 + *Ylik^k L e k^+' so ^ na ^ eacri atom is left in the excited state with a small probability. To 
initialize the state for the next trial, we allow the atoms to decay for a time T init . The total 
initialization error can then be estimated conservatively as the error in a single trial, 
exp(— Ti n it /r), times the number of trials. While the initialization step is by far the slowest 
step in an individual trial, the exponential scaling of n- m i t implies that the initialization 
time per trial scales only as logiV for entanglement of N particles. 



Implementations: neutral atoms and NV centres 

As an example, we consider the application of our method to a string of neu- 
tral atoms placed at adjacent antinodes of an optical lattice. This situation can be 
realized experimentally in the near future, due to recent advances in micromanipula- 



tion of neutral atoms 
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23]. Those experiments used Cs atoms in a lattice with 



532 nm spacing between antinodes. For circularly polarized laser excitation on the 
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ji\r = A,mp = 4) — > 2 P 3 / 2 |F = 5,mp = 5) cycling transition at 852 nm, the atoms 
behave as two-level systems. Laser excitation pulses of femtosecond duration and energies 
of a few nJ can implement essentially instantaneous n/2 rotations, and are readily obtained 
from mode-locked Ti: sapphire lasers. Temporally gated detection of the atomic emission at 



the single-photon level can be accomplished by upconverting the emission wit 
laser pulse and spectral filtering to detect only the up converted signal 
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i an auxiliary 
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Fig. [2] shows the emission pattern predicted by Eq. ([3]) for a string of 30 atoms lying 
parallel to The usual single-atom intensity distribution under <7 + -polarized excitation 
is Io(k) oc (1 + cos 2 #), where 9 is the angle between k and kp. Table [I] gives examples of 
entangled state preparation with 10, 30, and 100 atoms, with detection parameters chosen 
so that each error source contributes an error of less than 0.2/iV, ensuring the final state 
is entangled. We define S as the spontaneous emission rate in units of 1/r. Since S is 
approximately 1 in all cases, the atoms are not in the superradiant regime. The preparation 
times of a few ms are on the order of the time required for generation of six- and eight-ion 
entanglement in recent experiments 2lJ, [271]. Notably, the sum of the computed errors in 
preparing the 100-atom W state is considerably below the current gate error threshold for 
fault-tolerant quantum computation 281 ]. 



Our entanglement preparation method is useful for a broad range of qubit implemen- 
tations, in particular for ensembles of nitrogen-vacancy (NV) centres in diamond. The 
NV-centre defect exhibits a simple level structure with a zero-phonon line near 637 nm, and 
has been intensively studied for QIP applications 0, [|o]]. Optical pumping experiments 
with NV-centre ensembles have achieved up to 80% spin polarization of the ground 3 A 
state 3JJ, and inhomogeneous linewidths as low as 20 GHz have been observed on the 



zero-phonon line [32J]. The Amj = selection rule on the zero-phonon line ensures that 
initially sp in-p olarized centres behave as near-perfect two-level systems under optical 
excitation 33J], with dipole emission pattern Io(k) oc cos 2 9 for 7r-polarized excitation. 
Table [I] gives parameters for cylindrical ensembles of NV centres, using the intensity 
distribution derived for a cylinder by {34] and a value of 13 ns for the single-emitter decay 
time [33J. The volume of the cylinder is chosen to maintain a constant number density of 
2 x 10 14 cm -3 , and the length is taken to be ten times the diameter, yielding superradiant 
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TABLE I: Parameters for entangled state preparation with Cs atoms and NV centres, r, upper- 
state lifetime. N, number of qubits. S, spontaneous emission rate per atom in units of 1/r. a^, 
half-angle of detection cone, rydetj photon fraction emitted into detection cone. T^et, detection time. 
Ti n it, initialization time. N tT , number of trials to prepare entangled state with 50% probability. 
2p rcp , time required to prepare entangled state with 50% probability. 
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s, r = 


30 ns 


NV, r = 


13 ns 


N 


10 


30 


100 


1 n 


30 


100 


S 


1.2 


1.2 


1.2 


1.9 


3.0 


5.7 


adct [mrad] 


210 


91 


37 


280 


150 


62 


Vdet 


0.074 


0.040 


0.021 


0.17 


0.082 


0.026 


T det [ps] 


51 


5.6 


0.50 


1.6 


0.96 


0.046 


Timt [ns] 


255 


330 


400 


150 


160 


200 


iVtr 


470 


2600 


1.6 x 10 4 


1800 


1300 


1.3 x 10 4 




120 


840 


6400 


270 


200 


2800 



behavior (S ^> 1). As before, all detection parameters are chosen so that each error source 
contributes an error less than 0.2/N, ensuring the final state is entangled. Unlike the case 
of Cs atoms, the inhomogeneous broadening induces phase rotations between NV centres, 
but the accumulated phase is small for sufficiently short detection times. The resulting 
error is smaller than the other error sources for the examples in Table [H and limits the 
detection time only in the 10-qubit example. 



S 




FIG. 2: Angular intensity distribution of emitted radiation (red) for thirty Cs atoms in an optical 
lattice. The atoms (not shown), at the center of the plot, are equally spaced by 532 nm along 
the direction ki (green) of the exciting laser pulse. The intensity distribution is symmetric under 
rotation about k^. The collective forward emission peak is absent in the angular distribution for a 
single atom (black, not to scale). 



Conclusion 



We have presented a repeat-until-success protocol for entangling noninteracting qubits 
by detecting collective spontaneous emission. The protocol is applicable to a wide range of 
qubit implementations that exhibit spontaneous emission. The experimental setup is simple 
and robust, involving only laser pulses and single-photon detectors, with no need for optical 
resonators. This simplicity should enable practical realisation of mesoscopic entanglement 
in the near future. Built-in entanglement purification mechanisms let us reduce important 
error sources to any desired level by modifying the detection setup. Entangled states 
of at least 100 qubits can be generated with reasonable experimental parameters and 
with preparation times comparable to those required for generation of few-ion entangled 



states 



21 



271 ]. The protocol is especially useful for solid-state qubit implementations 



such as NV centers and quantum dots because of its resistance to inhomogeneous broadening. 
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